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Correlated transport of FQHE quasiparticles in a double- antidot system 
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Abstract 

We have calculated the linear conductance associated with tunneling of individual quasiparticles of primary quantum Hall liquids 
with filling factors v = l/(2m + 1) through a system of two antidots in series. On-site Coulomb interaction simulates the Fermi 
exclusion and makes the quasiparticle dynamics similar to that of tunneling electrons. The liquid edges serve as the quasiparticle 
reservoirs, and also create the dissipation mechanism for tunneling between the antidots. In the regime of strong dissipation, the 
conductance should exhibit resonant peaks of unusual form and a width proportional to the quasiparticle interaction energy U. 
In the weakly-damped regime, the shape of the resonant conductance peaks reflects coherent tunnel coupling of the antidots. The 
Luttinger- liquid singularity in the rates of quasiparticle tunneling to/from the liquid edges manifests itself as an additional weak 
resonant structure in the conductance curves. 
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1. Introduction 

Quasiparticles of two-dimensional (2D) electron liquids 
in the regime of the Fraction Quantum Hall effect (FQHE) 
have the unusual properties of fractional charge ^1^ and frac- 
tional exchange statistics |2l3j . Quantum antidots formed 
in the 2D electron system offer a possibility of localizing and 
controlling transport of individual quasiparticles [4|5|6j . 
Such a control made possible the first direct observation 
of the fractional quasiparticle charge in tunneling through 
an antidot ^J. This observation was later extended to the 
tegime of ballistic quasiparticle transport [718] . 

The mechanism of quasiparticle localization on antidots 
relies on the combined action of the electric and magnetic 
field and is microscopically quite different from the corre- 
sponding features of electron localization in quantum dots. 
Nevertheless, transport phenomena in antidots are very 
similar to those associated with the Coulomb blockade in 
tunneling of individual electrons 9 . For instance, in close 
analogy to the Coulomb-blockade oscillations of conduc- 
tance of the quantum dots [10 , the antidots exhibit peri- 
odic conductance oscillations with each period correspond- 
ing to addition of one quasiparticle to the antidot ^4^6] . So 
far, antidot transport has been studied both experimentally 
|4|5|6j and theoretically [TT| for one antidot. In this work, we 
develop a theory of correlated transport of individual quasi- 
particles through two antidots. The double- antidot system 



was discussed previously [121 as a qubit, information in 
which is encoded by individual quasiparticles. Such qubit is 
similar to superconducting charge qubits [13114] which are 
based on the dynamics of individual Cooper pairs. As in the 
case of Cooper pairs [15ll6j , the transport measurements on 
the quasiparticle qubit can be done more easily that direct 
measurements of the qubit dynamics. Transport measure- 
ments would constitute the first step towards experimental 
development of the FQHE qubits. More generally, under- 
standing the transport properties of multi-antidot systems, 
in particular the role of Coulomb interaction for localization 
of individual quasiparticles, and significance of the edge- 
state decoherence, should also be important for other, more 
complicated types of suggested FQHE qubits [17118] which 
also require control over individual quasiparticles. 

2. Model 

The system we consider consists of two antidots in se- 
ries between the two opposite edges of a primary quantum 
Hall liquid with the filling factor i/ = l/{2m+l) (Fig.[lji). 
The antidots are tunnel-coupled to each other and to the 
edges, which play the role of quasiparticle reservoirs. The 
quasiparticle current through the antidots is driven by the 
transport voltage V applied between the edges. The focus 
of this work is on the regime when all relevant energies are 
smaller than the energy gap A* of the antidots (see below) , 
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Fig. 1. Quasiparticle tunneling in the double-antidot system: (a) the 
real-space geometry (not to scale) of quasiparticle transfer between 
the opposite edges of the FQHE liquid; (b) energy diagram of the 
transfer. 

and the transport can be described completely in terms of 
the transfer of individual quasiparticles. This regime is rel- 
evant, e.g., for the operation of this system as a qubit. The 
main elements of the model of the double-antidot system 
in this case can be outlined as follows. 

2.1. Antidots 

An antidot formed at a point <^ in a primary quantum Hall 
liquid with the filling factor i' — 1/ {2m+l) can be described 
as a collection of n quasihole excitations created at this 
point. Microscopically, the unnormalized wavefunction of 
this configuration is [IJ: 

^{{^^A)-l[i^^,~Cr^m{{Zj}), (1) 

j 

where iprn{{zj}) is the Laughlin's wavefunction of the 
unperturbed quantum Hall liquid and, in the standard 
notations, the antidot position C in the two-dimensional 
plane and the electron coordinates Zj are given in the 
complex form and normalized to the characteristic length 
i = (n/eB)i/2 in the magnetic field B. The number n of 
the quasiholes is related to the geometric radius R of the 
antidot: R ~ ^/2n£. 

In what follows, we make use only of the general qualita- 
tive features of the wavefunction ([1]). For instance, in agree- 
ment with the typical experimental situation (see, e.g., [J), 
we assume that the antidot is relatively large: n 3> 1. This 
means that addition/removal of individual quasiparticles 
(here and below, this term will be used to describe pro- 
cesses with varying n: n — > n ± 1), does not change the 
antidot parameters noticeably. Indeed, the variation of the 
antidot radius in this case is 6R oc £/y^ and is small not 
only on the scale of R, but, more importantly, on the scale 
of the magnetic length £. 

The general form of the antidot energy £"„ as a function 
of n is determined by the interplay of Coulomb interaction 
and an external potential used to create the antidot. How- 
ever, for large n, and in some small range of variation of n 
around the minimum of En, one can always approximate 
this dependence as quadratic. This defines the character- 



istic energy gap A* = d^En/dv? which gives the energy 
interval of variation of the chemical potential /i of the sys- 
tem between the successive additions of individual quasi- 
particles to the antidot. For the system shown in Fig. [T^, 
the antidots exchange quasiparticles with the edges, and 

is defined by the edge chemical potential. In the situa- 
tion of the antidot, when all the energies are dominated 
by the Coulomb repulsion, the energy gap A* for chang- 
ing the number of quasiparticles is approximately related 
to the energy gap for the antidot excitations at fixed n: 
A* ~ hu/2'KR, where u is the velocity of the excitations en- 
circling the antidot. In general, e.g. in the case of quantum 
dots, the two types of energy gaps can be very different. 

We assume that the gap A* is sufficiently large for both 
antidots of the double-antidot system, so that in the rea- 
sonably large range of variation of both antidots are char- 
acterized by some well-defined numbers rt;, / — 1,2, of the 
quasiparticles. In this regime, the non-vanishing conduc- 
tance of the double-antidot system requires that /i is close 
to resonances in the both antidotes. At resonance, ~ 
Eni+i, and each antidot can in principle be in one of two 
states which differ by the presence/absence of one "extra" 
quasiparticle. The resulting four double-antidot states are 
relevant for the quasiparticle transport. We will use the no- 
tation for these states that gives the number of extra quasi- 
particles on each antidot: 

\i3) = \ni+i,n2+ i) , z,j = 0, 1, (2) 

and talk about the "first" and the "second" quasiparticle on 
the antidots disregarding the background ni quasiparticles. 
Counting the antidot energies En from the energy of the 
state with no extra quasiparticles, we can parameterize the 
energies tij of the four states as 

eoo = 0, eii = 2e + t/, eio = e - (5 , e^i ^ e + 5 . (3) 

Here U is the interaction energy between the extra quasi- 
particles on the two antidots, and e and 5 give the energies 

(Fig.[T)D) of the single-quasiparticle states localized at the 
two antidots: £1^2 = e ± (5. The energies e and 5 are defined 
relative to the common chemical potential of the edges for 
vanishing bias voltage V between them. Non- vanishing bias 
voltage shifts the energies ([3]) . Experimentally, the antidot 
energies are controlled by the back-gate voltage or mag- 
netic field [415] . The degree to which these fields couple to 
the energy difference 5 depends on the degree of asymme- 
try between the two antidots. In the following, we present 
the results for quasiparticle conductance of the system as 
a function of e for fixed (5, as would be appropriate for 
identical antidots. These results can be generalized to non- 
identical antidots by taking a "cross-section" in the space 
of e and 5 along the direction appropriate for a given degree 
of the antidot asymmetry. 

If the two antidots are sufficiently close, so that the dis- 
tance between their edges is on the order of magnetic length 
I, the quasiparticle states localized around them overlap 
and hybridize. This effect can be accounted for by the tun- 
nel coupling — A of the antidots. The phases of the antidot 



states can always be chosen to make A real. This coupling 
affects only the singly-occupied states |10) and |01). The 
single-quasiparticle part of the Hamiltonian is then: 



H 



5(7^ — AcTt 



(4) 



where cr's are the Pauli matrices. Equation (jl]), together 
with the part of Eq. ^ describing states with zero and 
two quasiparticles, give the main part of the antidot energy 
controlling the quasiparticle transport. In what follows, we 
assume that all contributions to this energy and the tem- 
perature T are small: 



A,e,,[/,T« A* 



(5) 



In this regime, tunneling through the antidots can be dis- 
cussed in terms of correlated transfer of individual quasi- 
particles. Since the gap A* is dominated by the Coulomb 
interaction, i.e. has the same origin as the quasiparticle in- 
teraction energy U , the most restrictive part of the assump- 
tion ([S]) is the condition on U . This condition can still be 
satisfied, due to the difference between stronger quasipar- 
ticle repulsion on the same site and weaker repulsion on 
different sites. 

2.2. Antidot- edge tunneling 

Similarly to the tunnel coupling between the antidots, 
if the edges of the FQHE liquid are not far from the anti- 
dots on the scale of the magnetic length i, there is a non- 
vanishing amplitude for quasiparticle tunneling between 
the edge and the nearest antidot. The tunneling between 
the lt\i edge and antidot can be described quantitatively 
with the standard tunnel Hamiltonian 



Hr 



Titpki + h.c. 



(6) 



where '0,'0^ ^ind ^i^^ &re the creation/ ahhihalation oper- 
ators for quasiparticles at, respectively, the edges and the 
antidots. Denoting the position along the edge as x and 
taking the tunneling points for both edges to be at x = 0, 
the edge quasiparticle operators ijji can be expressed in the 
standard bosonisation approach as [19] 



4>i(t) = (l/27ra)i/2|^g.v^0,(o,t) 



(7) 



Here the "Klein factors" ^; account for the mutual statis- 
tics of the quasiparticles in different edges, 0/ are the chiral 
bosonic fields which describe the edge fiuctuations propa- 
gating with velocity u, and 1/a is their momentum cut-off. 
The edge fluctuations result in the fluctuations of electron 
density at the edge: pi{x,t) — {^/v /2T:)d(f)i{x,t) / dx. Both 
fields (pi can be decomposed in the standard way into the 
individual "magneto-plasmon" oscillator modes a„,ajj: 



n— 1 



, g„ = 27rn/L , (8) 



where L is a normalization length. In the geometry of the 
double-antidot system, there are no interference loops for 



the edge quasiparticle. In this case, the statistical Klein fac- 
tors in ^ would cancel out in the perturbation expansion 
in tunneling ([6]) and can be omitted. 

The quasiparticles at the antidots should be described 
in general by the expressions similar to Eq. ([7]). Condition 
([5]) of the large antidot energy gap A* ensures, however, 
that the fluctuations of the edges around the antidots are 
suppressed, i.e. the magneto-plasmon oscillations are not 
excited out of their ground state |0). In this regime of the 
"quantized" edge, the general quasiparticle operators ([7]) 
reduce to just the statistical Klein factors up to a normaliza- 
tion constant. Indeed, as one can see directly from Eq. ([5]) 
by bringing the (/)-part of ^ into the normal form. 



s''^^|0)/(2™)i/2 = (7ri?)-i/2 



(9) 



Including this normalization constant in the tunnel ampli- 
tude T/, we see that the operators ^; for quasiparticles at 
the antidots consist solely of the Klein factors. The appro- 
priate set of properties of the quasiparticle Klein factors ^ 
depends on the specific geometry of each edge-state tun- 
neling problem. Non-trivial examples of this can be found 
in |20|21|22|23l24] . As follows from the discussion in the 
preceding Section, in tunneling between the quantum an- 
tidots, the operators should account for the "hardcore" 
property of the quasiparticles: In the given range of exter- 
nal parameters only one extra quasiparticle can occupy one 
antidot. In general, these operators should also describe 
the anyonic exchange statistics of the FQHE quasiparticles, 
but the geometry of the double-antidot system (Fig. [T^) 
does not permit quasiparticle exchanges, and the exchange 
statistics of the tunneling particles is irrelevant. Since the 
hardcore property makes the quasiparticle occupation fac- 
tors equivalent to those of the fermions, and the actual 
exchange statistics is irrelevant, the antidot quasiparticle 
operators ^, can be treated as fermions. Together with 
Eqs. ([7]) and ^ for the edge quasiparticles, this defines 
completely the tunnel Hamiltonian ([6]) . 



2.3. Edge- state decoherence 

Tunneling of charged quasiparticles through the antidot 
system couples to all gapless charged excitations that ex- 
ist in the system. In the case of the FQHE liquid, excita- 
tions in the bulk of the liquid are suppressed by the en- 
ergy gap, and only the edges support gapless excitations. 
In contrast to all other possible mechanisms of decoherence 
(e.g., plasmons in metallic gates, or charged impurities in 
the substrate) the edges play the role of reservoirs in trans- 
port measurements and as a matter of principle can not 
be removed from the antidots. In this Section, we estimate 
the strength of this unavoidable edge-state decoherence for 
quasiparticle tunneling through the double-antidot system. 

The spectrum of the gapless edge excitations of one 
edge consists of magneto-plasmon oscillations ^ with the 
Hamiltonian: 



Q 



Ho = [hu/L) J2 



(10) 



We assume that the antidot system is symmetric, and a 
quasiparticle sitting on the first antidot creates a potential 
Vi{x) along the lih edge. The quasiparticle dynamics gov- 
erned by the Hamiltonian ^ is coupled then to the fluctu- 
ations of electron densities pi (x) at the edges through the 
interaction Hamiltonian 



Hint = ^cr^e 



dxV{x){pi{x) - P2{x)) , 



(11) 



where V{x) = Vi{x) — V2{x) is the change of the potential 
along the edges due to quasiparticle transfer between the 
antidots. Since the edge-antidot distance and the distance 
between the antidots are on the order of antidot radius i?, 
this radius sets the range of the potential V{x). The edge 
velocity u can be expected to be similar for the external 
edges and the antidots. This means that the condition ([5]) 
of the large energy gap implies that the characteristic wave- 
length of the edge excitations which can exchange energy 
with the quasiparticles on the antidots is much larger than 
the range of the potential V{x): hu/e ^ hu/A* ~ R. The 
interaction energy Ijlip can be expressed then as 



H^nt = ^<yz{pi{0) - P2(0)) j dxVix) , 

where, as follows from Eq. ^ , the densities pi are 

. , — C30 

P;(0) = ^ Vn{an - 4) • 

n=l 



(12) 



(13) 



The strength of interaction (|12p can be characterized by 
the typical transition rate between the eigenstates of the 
antidot Hamiltonian ([4]) induced by the edges. Straightfor- 
ward calculation of the "Golden-rule" rate using Eqs. ([TT 
(HD, and (flHl gives: 



AE 



a K {a. 



1 



-AE/T 



(14) 



where AE is the energy difference between the two states, 
(az) is the matrix element of CTz between them. The dimen- 
sionless factor k characterizes the overall "intensity" of the 
antidot-edge potential: 



yAneeoJ 



dxV{x) . 



(15) 



The precise form of the potential V{x) and the value of 
K, depend on the details of configuration of the metallic 
gates that define the edges and screen the antidot-edge 
interaction. However, normalized as in Eq. (jlSp . n should 
be on the order of 1. For instance, assuming as a crude 
model of the system electrostatics that the antidot-edge 
interaction is confined to the interval d ~ 2R in which the 
edge is a tunnel- limited distance £ away from the antidot, 
one can estimate k as {2/tt) ln(2R/£), i.e., k ~ 2 for realistic 
10. 



The factor a in (|14p is the "fine structure constant" of 
the edge excitations: 



(16) 



Aneenhu 



and is the main parameter controlling the strength of deco- 
herence Td through the velocity u of the edge excitations. 
The dielectric constant e is fixed by the material (GaAs) 
of the structure: e ~ 10, and in the realistic range of pos- 
sible velocities u, 10'* 10^ m/s [S], a should vary in the 
range between 2 and 20. In the most relevant case of the 
FQHE liquid with the filling factor ly = 1/3, and for the 
edge-antidot coupling intensity estimated above, the qual- 
ity factor AE/hTd of the quasiparticle dynamics changes 
then roughly between 0.1 and 10. This means that in the 
case of strong edge confinement that produces large ve- 
locity u, the quasiparticle dynamics on the antidots can 
be quantum-coherent provided that all other decoherence 
mechanisms are sufficiently weak. In the opposite case of 
smooth confinement with low velocity m, the already un- 
avoidable edge-state decoherence is strong enough to com- 
pletely suppress the coherence of the quasiparticle states 
on different antidots, and quasiparticle transfer processes 
between them are incoherent. 



3. Tunneling rates 

As was mentioned above, the discussion in this work is 
limited to the regime in which the transport through the 
double-antidot system can be interpreted as the correlated 
transfer of individual quasiparticles. Besides the condition 
(O on antidot energies, this also requires that the antidots 
arc coupled only weakly to the edges, so that the edge- 
antidot tunneling can be treated as a perturbation leading 
to an incoherent transfer of individual quasiparticles. The 
quasiparticle transport through the antidots is governed 
then by the kinetic equation similar to that for Coulomb- 
blockade transport in quantum dots with discrete energy 
spectrum ,25J . This Section calculates the relevant tunnel- 
ing rates in the two limits of strong and weak edge-state 
decoherence. 



3.1. Strong decoherence 

If the edge-state decoherence is sufficiently strong, the 
quasiparticle transfer between the antidots can be treated 
as incoherent and described by the sequential tunneling rate 
obtained by perturbation theory in the tunnel amplitude A. 
To calculate this rate, it is convenient to express the density 
operators pi of the two edges through one effective density 
p which satisfies the same relation ([T3|: pi{Q) — /O2(0) 
-\/2p(0), so that the edge-antidot coupling (fT2|) is: 



Hint = ^^^^ P^^"^ 



(17) 



A 



Next, one can perform a unitary transformation which con- 
verts the fluctuations of the energy of the quasiparticle ba- 
sis states 1 10) , |01) induced by p7|) into a fluctuating phase 
of the tunneUng matrix elements of the tunneling part of 
the Hamiltonian (jl]): 



27r2 



(18) 



where ^(x, t) is the bosonic field given by same the Eq. dS]). 
Then, the rate Fa of the anidot-antidot tunneling can be 
expressed in the lowest non-vanishing order in the ampli- 
tude A as: 



2A2Re / dte 



iEt 



/giv^0(O,t)g- 



(19) 



where (...) is the average over the equilibrium fluctuations 
of (/) and E — ±25 is the energy difference (depending on the 
direction of tunneling) between the quasiparticle states lo- 
calized on the antidots. The standard evaluation of Eq. 
(see, e.g., EH?]) gives: 



T^{E) = -ifg{E), 7 = 27rA2K 



2TiT{g) 



{2t:T/uj^Y-^ \T{g/2 + iE/2TiT)\^ e 



(20) 



-E/2T 



where T{z) is the gamma-function and lOc = hu/2a is 
the cut-off energy of the edge excitations. The function 
fg{E) gives the energy dependence of the tunneling rate 
(see Fig. [5]) and is defined to coincide with the Fermi dis- 
tribution function for g = 1. The power g determines the 
behavior of the transition rate at large energies \E\ ^ T: 
Fa(£^) oc E^^~^^ on the "allowed" side of the transition 
{E < 0), and Fa(£^) cx E^s-i)^-!^/^ on the "forbidden" 
side {E > 0), when the transition has to overcome the en- 
ergy barrier E. This asymptotic behavior of the tunneling 
rates, together with Eq. ((^ . is valid at \E\ <^ lOc- 

The rates F;, I = 1,2, of the antidot-edge tunneling are 
obtained through a similar calculation starting with the 
tunnel Hamiltonian ([5]). They are given by the same ex- 
pression (pO| : 



TiiE)=jiME), ji 



2^|TipK 



(21) 



In general, the long-range Coulomb interaction should gen- 
erate corrections to g which move it away from the "quan- 
tized" value g = v [28j . However, in contrast to the quasi- 
particle tunneling between the antidots, which is changed 
qualitatively by decoherence created by the Coulomb in- 
teraction with the edge, the Coulomb corrections for the 
antidot-edge tunneling are expected to be small and will 
be neglected in this work. 

Thus, in the regime of strong edge-state decoherence, 
the overall transport of quasiparticles through the double- 
antidot system can be described as a combination of suc- 
cessive antidot-edge transitions (|2T]) and incoherent tran- 
sitions (PO)) between the antidots. 




Fig. 2. Energy dependence of the antidot-antidot II20I I and an- 
tidot-edge II21I I tunneling rates. The normalization factor is 
To =7(27rT/a;c)9-i. 

3.2. Weak decoherence 



For sufficiently strong edge-state confinement, the edge- 
induced relaxation rate (fT4|) will be smaller than the antidot 
energies. If other decoherence mechanisms, including the 
decoherence created by incoherent antidot-edge tunneling 
are also weak on the scale of the antidot energies. 



Fd,Fi « A,,S,f/,T, 



(22) 



the quasiparticle dynamics on the antidots is quantum- 
coherent. It is characterized by the stationary eigenstates 
\k), k = 1, 2, of the double-antidot Hamiltonian (|4]): 

H\k) = {e + {-l)''n)\k), n={S^ + A2)1/2^ 

|fc) =cifc|10)+C2fc|01), (23) 

for which the probabilities Xik of finding the quasiparticle 
on the Zth antidot are: 



Xik^\cik\^ = [i + {-iy+'s/n]/2 



(24) 

In the coherent regime (P^ . the double-antidot system 
can be viewed as a quasiparticle qubit [H]. The current 
through the qubit is described in terms of tunneling to/ from 
the eigenstates ([23l) . The corresponding tunneling rates are 
found from the tunnel Hamiltonian ([6|) , in which, as was dis- 
cussed in Sec. I2.2| the quasiparticle creation/annihalation 
operators ^, act as fermions. This means that the tunnel 
matrix elements for the quasiparticles can be calculated in 
the standard way. In particular, for each eigenstate ([23]) . 
the matrix element is independent of the occupation factor 
of the other eigenstate. Explicitly, the tunneling rate Tik 
from the Ith edge into the state |fc) is 



Ik 



2|T,n(fc|^/|0)pRe / dte^^\^j{t)MO)) 



(25) 



Here |0) denotes the empty eigenstate and E is the ap- 
propriate tunneling energy which includes in general the 



eigenenergies (j23|) and the interaction energy U. The quasi- 
particle matrix elements are |(fc|f/|0)p — Xik, and we get: 

TikiE) ^ XikTiiE), (26) 

where the rates Ti{E) are given by Eq. ([2T|l . In the prac- 
tically important case of FQHE liquid with v — 1/3, the 
energy dependence of the transition rates ([26| is illustrated 
by the g — 1/3 curve in Fig. [51 The peak of the tunnel- 
ing rate at e ~ is the consequence of the Luttinger-liquid 
correlations of the edge quasiparticles. Conductance calcu- 
lations presented in the next Section show that this peak 
manifests itself as additional resonant features of the qubit 
conductance. 

4. Conductance of the double-antidot system 

In both situations of strong and weak decoherence, the 
conductance associated with tunneling of individual quasi- 
particles through the double-antidot system can be calcu- 
lated by solving the kinetic equation for quasiparticle occu- 
pation probabilities of the antidot states. Similarly to the 
case of tunneling through one antidot |4|5|6|llj , the con- 
ductance as a function of the common energy e of the anti- 
dot states should exhibit the resonant tunneling peaks. For 
the double-antidot system, the peak structure is, however, 
more complicated, reflecting the transition between the 
low-temperature regime in which each peak corresponds to 
addition of one quasiparticle to the system of antidots, and 
a possible "large-temperature" regime, when the single- 
quasiparticle peaks are merged, and each conductance peak 
is associated with addition of two quasiparticles. In this 
Section, we calculate the corresponding conductance line 
shapes. Quantitatively, these line shapes are determined by 
the interplay between the quasiparticle repulsion energy U 
on the two antidots and tunnel coupling A between them. 
The calculations below are focused mostly on the more typ- 
ical case of large repulsion energy U ^ A ~ S. 

4.1. Strong decoherence 

For strong edge-state decoherence, coherent mixing of 
the quasiparticle states on the two antidots is suppressed, 
and the quasiparticle dynamics is described by kinetic equa- 
tions for the occupation probabilities pij of the states \ij) 
^ localized on the antidots. The probabilities evolve due 
to incoherent jumps of quasiparticles at the rates Fa (I20| 
and F/ ((2T|) between these states. The stationary quasipar- 
ticle current / through the antidots is found in this regime 
from the balance of the forward/backward transition across 
any of the three tunnel junctions of the system, e.g., from 
the transitions between the antidots: 

/ = ez/[pioFA(2(5) - poirA(-2(5)] . (27) 

In general, the quasiparticle current / can be calculated 
by the direct numerical solution of the kinetic equation. 



The results of such solution for the linear conductance G = 
dl / dV\v=Q are shown in Fig. [3] (in all numerical results 
presented below we take v = 1/3). Qualitative behavior of 
the system can be understood by analyzing the limits where 
the simple analytical expressions for the conductance can 
be obtained. 

The first limit is Fa ^ F; , where the antidot-antidot tun- 
neling is the bottleneck for the current flow. In this case, to 
the zeroth-order approximation in Fa, the current is van- 
ishing, and one can use in Eq. (I27|) the equilibrium proba- 
bilities = (1/Z)e~'^'5/'^, Z — J2ij e~''j/^, obtaining for 
the conductance: 

T e-^/^ + e-(^+t^)/^ + 2cosh(VT) ■ ^ ' 

Equation ([28]) describes the "coalesced" conductance 
peak that corresponds to the addition of two quasiparti- 
cles to the antidots. At large temperatures, T > U, the 
peak has a usual thermally-broadened shape with width 
proportional to T. At T <C U, however, the peak shape 
(|28p is quite unusual: the conductance is constant between 
the point e ~ 0, when the first quasiparticle is added to the 
antidots, and the point e ~ —U of addition of the second 
quasiparticle, forming the plateau of width U - see Fig. [31 
The conductance plateau remains flat until the tempera- 
ture is lowered to T ~ U / \n[Ti{U) /T a], when the thermal 
suppression of the antdot-edge tunneling rate makes it 
comparable to Fa at the center of the plateau, e ~ ~U/2, 
despite the fact that the two rates are very different at 
e ~ 0. In this temperature range, a dip develops in the 
center, which separates the plateau into two peaks, one at 
e ~ and the other at e ~ — ?7, with decreasing temper- 
ature (Fig. [3]) . Each peak corresponds to addition of one 
quasiparticle to the double-antidot system. Note that the 
resonant peaks occur when the gate bias energy e is equal 
to "minus energy" of the antidot state, so that the total 
energy of the state relative to the chemical potential of the 
edges is zero. 

The shape of such single-quasiparticle peaks can be de- 
scribed in the opposite limit of strong antidot-antidot tun- 
neling Fa 3> Fj. In this limit, the general kinetic equations 
for three probabilities poo,pio,poi relevant at e ~ 0, e.g.. 

Poo = ri(-ei)pio + F2(-e2)poi - [ri(ei) + F2(e2)]poo, (29) 

and similar equations for the other probabilities, can be 
reduced to two equations for the effective two-state sys- 
tem. The strong antidot-antdot tunneling that couples 
the singly-occupied states |10),|01), maintains the rela- 
tive equilibrium between them: pw/poi = e^^*/^, making 
it possible to treat these two states as one. The effective 
transition rates between this state and the state |00) are 
obtained as the weighted average of the transition rates in 
starting kinetic equations, e.g. (|29|) . The standard calcula- 
tion of the current through a two-state system gives then 
the resonant peak of the double-antidot conductance at 
e ~ (associated with addition of the first quasiparticle to 
the antidots): 




Fig. 3. Conductance of the double-antidot system in the regime of 
overdamped quasiparticle transport. Conductance is normaUzed to 
Go = (ej^)^r/^(2(5)/T. Different curves correspond to different ratios 
of the antidot-antidot and antidot-edge tunnehng rates. From top to 
bottom: rA(2(5)/r2(0) = IQ-*, 10"^, 10"^, 10"^ 
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where f{E) = fi{E) is the Fermi distribution function. For 
instance, if the two antidot states are ahgned, S — 0, 
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The conductance peak at e ~ — J7 associated with the 
addition of the second quasiparticle is given by an expres- 
sion similar to Eq. (|30p with an appropriate shift of energy 
e ^ e + J7. In particular, for 5 = 0, this expression reduces 
to 

evf 71 72 f„{e + U) 



G 
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The conductance peak (PTjl is asymmetric around e = 0, 
since in the tunneling dynamics underlying this peak, only 
one quasiparticle can tunnel off the antidots, while there 
are two available states for tunneling onto the antidots. 
Still, the "quasiparticle-quasihole" symmetry makes the 
two peaks, (fSTj) and ([32|l . at S = symmetric images of 
each other with respect to a "mirror" reflection e + U/2 ^ 
-(e + U/2). For (5^0, the condition Fa > F/ is violated 
at sufficiently low temperatures T <C (5, and Eq. ((30|) be- 
comes invalid. In this case, the antidots are effectively out 
of resonance, and conductance peaks are suppressed expo- 
nentially with temperature at all gate bias energies e. 

4.2. Weak decoherence 

If the edge-state decoherence is sufficiently weak and al- 
lows for quantum-coherent transfer of quasiparticles be- 
tween the two antidots, the kinetic equation for quasipar- 
ticle transport should be written not in the basis of states, 
([2]), but in the basis of the hybridized states (|23|) . As fol- 
lows from the estimates of the edge-state decoherence in 



Sec. 12.31 even in this regime, the edge-induced relaxation 
rate Td (|14p should be strong enough, F^j ^ F;, to main- 
tain the equilibrium distribution of quasiparticles over the 
antidot states in the process of tunneling. This means that 
if i?^"-* is the energy of the state \k) when there are n quasi- 
particles on the antidots, the probability that this state is 
occupied is pkin) = (l/Z„)e-^i"'/^, Z„ = Efce"^^"'/^. 
The quasiparticle tunneling is reduced then to the dynam- 
ics of the total number n of quasiparticles on the antidots, 
described by the probability distribution p(n). The rates of 
tunneling transitions n ^ n ± 1 in this dynamics are: 

r^W= E^i^N' r±(n) = ^pfc(7i)F,(fc,g,n,n±l). 

1=1,2 kq 

where the partial transitions rates Ti{p, k,n,n ± 1) from 
the state k of n quasiparticles into the state A:' of n ± 1 
quasiparticles are given by the appropriate tunneling rates 
(|26|) between the ^th edge and antidot. The solution of the 
simple kinetic equation 

p{n) = [r^(n ± l)pin ± 1) - F±(n)p(n)] , (33) 



gives then the stationary quasiparticle current through the 
system: 



/ = jye ^ [F+ (n)P(n) - F^; (n - l)P{n - 1)] 



(34) 
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Fig. 4. Conductance of the double-antidot system in the regime of 
the underdamped quasiparticle dynamics. Conductance is plotted in 
units of Go = (ej/)-^ri(0)/T. The curves show the two main resonant 
conductance peaks at e = H and e = —{U + H), and a weak kink 
at e = —Q that is made visible by the Luttinger-liquid singularity 
in the tunneling rates. The upper and lower curves are, respectively, 
the conductance with and without equilibration on the antidots. 



Equation (j33p shows that the stationary probability dis- 
tribution p{n) satisfies the "detailed balance" condition 
p{n)T~^{n) = p{n + l)F~(n -|- 1) even in the presence of 
the non- vanishing bias voltage V. Using this condition, and 



expanding both p{n) and the tunnehng rates T{n) to first 
order in V, one finds the linear conductance G of the quasi- 
particle qubit: 



T 



'r+(n) + r+H 



(35) 



Here w„ = Zn/Z is the equihbrium probabihty to have n 
quasiparticles on the antidots, Z — Z„, and the factor 
T] gives the fraction of the voltage V that drops across the 
edge-antidot junctions. Equation (j35p can be understood in 
terms of forward jumps of quasiparticles in the left junction 
contributing to the current only if they are followed by the 
forward jumps in the right junction. As an example, at 
temperatures T <^ U, and e ~ —Q one can limit the sum 
in Eq. (|35p to one term n 0. The conductance G is then: 



G = 
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l + 2e-^/^cosh(17/r) 



E,,feri,(e+(-l)^»)r2fc(6 + (-l)fc») 
where the tunneling rates Tqk are defined in Eq. ((26 
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Fig. 5. Conductance of the symmetric [S = 0) antidot qubit exhibit- 
ing two resonant conductance peaks at e = A and e = — (t/ + A). 
Both peaks have kinks at e = — A and e = —{U — A) caused by the 
Luttinger-liquid singularity in the tunnehng rates. The inset shows 
the conductance for the special value of interaction energy U = 2 A, 
when the two kinks coincide producing very small but visible addi- 
tional conductance peak. Conductance is normalized as in Fig. |4] 

For comparison, one can calculate the conductance in the 
same regime T ^ [/, e ~ —17, but without equilibration 
on the antidots, i.e. assuming that the edge-state decoher- 
ence is very weak, <C F;. As before, the antidots can be 
occupied in this regime at most by one quasiparticle at a 
time, and straightforward solution of the kinetic equation 
describing the occupation of individual energy eigenstates 
due to transitions gives the conductance: 
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(37) 



T 2172 1 + 2e-^IT cosh(rj /T) 



Equation ||T7|) describes the resonant conductance peak 
that corresponds to the addition of the first quasiparti- 
cle to the antidot. The second quasiparticle peak at e = 
— ([/ -I- 17) is described by the similar expression. At low 
temperatures, T <C 17, only the lowest energy eigenstate 
with energy —17 contributes to the conductance ([37|) . In 
this case, the equilibration on the antidots does not have 
any effect, and Eqs. ([36]) and ((37|) coincide. As one can see 
from Fig. [4l which plots the conductance obtained by nu- 
merical solution of the full kinetic equation, the difference 
between the two regimes, with and without relaxation, re- 
mains very small even at moderate temperatures. At larger 
temperatures, 17 ^ T ^ ?7, and J = 0, Eq. ([57)1 reduces 
to Eq. ((3T|l for the conductance in the overdamped regime. 
The only difference between the two results is the factor 
77 in Eq. ([57|) which implies that the part of the applied 
bias voltage that drops across the region of the quantum- 
coherent quasiparticle dynamics does not contribute to the 
linear conductance. 

Besides the two main resonant peaks, the curves in Fig.[?] 
exhibit also a small kink at e ~ —17. This kink appears 
at the intermediate temperatures and is the result of the 
transfer of the first quasiparticle added to the antidots not 
through the more probable ground state of the qubit but 
through the excited state with energy 17. One could see, 
however, by plotting the conductance of the double- antidot 
system for tunneling electrons (the tunneling rates given 
by the g — l\n Fig. ^ that the contribution of the excited 
state to the conductance is not sufficient by itself to produce 
such a kink. The kink in the conductance appears only when 
the contribution from the excited state is amplified by the 
Luttinger-liquid singularity in the quasiparticle tunneling 
rate (seen in the g ~ 1/3 curve in Fig. [5] as a peak at 
zero energy). It becomes somewhat more pronounced in 
the conductance peaks of the "symmetric" qubit with 5 = 
shown in Fig. [S] In this case, the kinks appear on both 
peaks: at e = —A and e = —([/—A). The second kink is due 
to transport through the ground state of the qubit in the 
regime when the main contribution to conductance comes 
from the excited state. As shown in the inset in Fig. [5l at 
the special value of the interaction energy U ~ 2A, the 
two kinks coincide and form a weak additional peak of the 
qubit conductance. 

5. Conclusion 

We have calculated the linear conductance G of the 
double-antidot system in the regime of weak quasipar- 
ticle tunneling through the antidots. Depending on the 
strength of the edge-state decoherence, the tunneling can 
be coherent or incoherent. In the incoherent regime, the 
two resonant conductance peaks that correspond to the 
two antidot states are spaced by the quasiparticle inter- 
action energy U . In the coherent regime, this spacing is 



Q 



increased toU + 2f2, where 2^ is the gap between the en- 
ergy eigenstates of the double-antidot system. The coher- 
ent regime of quasiparticle dynamics is also characterized 
by the Lorentzian dependence of the system conductance, 
G oc (1 + on the energy difference S between 

the antidots. In the quantum-coherent regime, the double- 
antidot system can be used as a quasiparticle qubit. 
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